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ABSTRACT 

A model for the evolution of low-luminosity radio galaxies is presented. In the model, the lobes in- 
flated by low-power jets are assumed to expand in near pressure-balance against the external medium. 
Both cases of constant external pressure and decreasing external pressure are considered. Evolution 
of an individual source is described by the power-size track. The source appears as its lobe is in- 
flated and radio luminosity increases to above the detection level; the source then moves along the 
track and eventually disappears as its luminosity drops below the detection limit. The power-size 
tracks are calculated including the combined energy losses due to synchrotron radiation, adiabatic 
expansion, and inverse Compton scattering. It is shown that in general, the constant-pressure model 
predicts an excess number of luminous, small-size sources while underpredicting large-size sources in 
the power-size diagram. The predicted spectra are steep for most sources, which is inconsistent with 
observations. By comparison, the pressure-limiting model fits observations better. In this model, 
low-luminosity sources undergo substantial expansion losses in the initial phase and as a result, it 
predicts fewer luminous, small-size sources. The resultant spectra are flat for most sources except for 
the oldest ones, which seems consistent with observations. The power-size tracks, in contrast to that 
of high-luminosity radio galaxies, are characterized by a slow increase in luminosity for most of the 
source's life, followed by a rapid decline when the synchrotron or inverse Compton scattering losses 
set in. 

Subject headings: acceleration of particles-radiation mechanisms:nonthermal-galaxies:active 
-galaxies: jets 



1. INTRODUCTION 

The generic model for radio galaxies assumes that twin 
jets emanating from an active galactic nucleus propa- 
gate outward in two opposite directions. The jets, which 
initially propagate at a relativistic speed, interact with 
the surrounding medium leading to formation of a dif- 
fuse emission region. Radio galaxies appear to have two 
classes: low- and high-luminosity radio galaxies, com- 
monly referred to as F R I and II sources, respectively 
(Fanaroff & Rilev 1974). The jets in high-luminosity ra- 
dio galaxies have relatively homogeneous morphology; 
they are well collimated and propagate through the sur- 
rounding medium — initially in the cores, then halos of 
their parent galaxies and then the intergalactic medium 
(IGM) — creating pair of large lobes. The jets are dim un- 
til the end of the lobes where there are bright hot spots. 
Classical double radio sources are a typical example of 
this class. 

By contrast, low-luminosity radio galaxies are charac- 
terized by jets that are bright close to the nucleus of 
their parent galaxy. The jets have diverse morphologies, 
a feature that can be interpreted as deceleration of jets 
due to entrainment of the external medium. The jets are 
initially laminar near the nucleus and then subject to tur- 
bulent disruption when passing through the flare region 
that is thought to be the main acceleration site for rela- 
tivistic particles. The jets beyond the flare region spread 
out, resembling smoke arising from a chimney mixing 
with the ambient medium. 

The key issues in the understanding of radio galax- 
ies include the evolution of radio galaxies and the un- 
derlining physics that distinguishes these two classes. 
One suggestion is that these two classes of source are 



intrinsically different, primarily in their jet dynamics , 
evolving along different tracks ( Jackson fc Walll I1999T ). 
However, there are suggestions that some of the high- 
luminosity radio sources with weak jets may evolve into 
low-luminosity sour ces JGopal-Krishna fe Wiital H987L 
1988; Kaiser fc Best 2007). Some radio sources exhibit 
mixed features of FR Is and lis. For example, there are 
souces with one-side jet showing the FR I features and 
the other showing the FR II features. This leads to an 
opinion that such c lassification may not be clear cut as 
previously thought (K aiser fc Bestlf2007f) . 

It is well accepted that the radio emission in radio 
galaxies is due to synchrotron radiation by relativistic 
electrons (or positrons) injected from the jets. The total 
synchrotron power P v evolves with time as the injection 
of a mixture of kinetic energy and magnetic energy com- 
petes against the losses due to volume expansion and 
radiation. When the losses dominate, the total power is 
a decreasing function of time as the source ages. Since 
the typical evolutionary time scale is ~ I0 8 yr, it is not 
practical to measure how the total power changes in time 
directly by observations. One may study the temporal 
evolution of radio galaxies from the total spec tral power, 
P v , a s a function of the source's linear size ()Shklovskiil 
196^). The linear size here is defined as the dimension 
of the lobe along the jet axis. Since the linear size D 
increases as the source expands, a radio source should 
evolve along a particular track in the P v -D diagram. 

There are many discussions in the literature 
on the time evolution of hig h-luminosity radio 
galaxies (or FR II sour c es) (iKaiser et al.l [T997t 
Blundell. Rawlings fc Willott] fl999l iManolakou fc Kirkl 



2002h . In the existing models, there are three rele- 
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vant regions where the physical processes determine the 
evolution of the source. These include the hot spots 
where particles are assumed to be accelerated and ra- 
diate, the head region that contains the hot spots, and 
the lobe — an emission volume inflated by the input of 
the jet, where relativistic particles are injected and the 
volume of the emitting plasma expands. Since the hot 
spots appear at the end of the jet, their distance to the 
center can be identified as the linear size of the lobe. 
Thus, the modeling of the time evolution of the lobe 
size is reduced to the problem of modeling of changes 
in the location of the hot spots. One of the widely- 
discussed models is the self-similar expansion model in 
which the jet cr eates a bow shock by sweeping up the am- 
bient material (IFedorenko fc Zentsovalll986l iFalld 119911: 
iKomissarov fe Falldll998l ). So. the location of the shock 
is completely determined by the jet power and the den- 
sity of the surrounding medium. This allows one to es- 
tablish the source size as a function of time. 

In this paper we consider the evolution of low- 
luminosity radio galaxies (or FR Is). We derive the radio 
power as a function of the source's size, P U (D), which 
can be directly compared with the power-size diagram 
inferred from observations. In contrast to FRIIs, there 
are few discussions on the evolution of the low-luminosity 
sources, mainly because of the lack of quantitative mod- 
els that relate the linear sizes to the jet dynamics. One 
of the defining features of FR Is is that the flare region, 
the hot-spot equivalent as compared to high-luminosity 
sources, is located close to the nucleus, which indicates 
that the jet is decelerated to the subsonic flow regime 
relatively close by the nucleus. Since the jet continues 
to expand well beyond the flare region, forming a diffuse 
emission region further beyond, the size of a FR I source 
is not directly related to the location of the flare region 
and should be determined separately from the expansion 
of the diffuse region. Here we continue to refer to this dif- 
fuse emission region as the lobe as in FR lis despite the 
significant difference between their morphologies. For 
low-power jets, the physical conditions of the surround- 
ing medium play a critical role in defining the change in 
the lobe size. If the medium is warm, the expansion pro- 
ceeds at near pressure balance against the ambient pres- 
sure and the volume increases more slowly with time than 
the self-similar expansion in the high-luminosity sources. 
As a result, the low-luminosity sources grow in size, by 
comparison, much more slowly than the high-luminosity 
sources. 

In Sec 2, we discuss a generic model for both high- 
power and low-power jets. The evolution of the emitting 
plasma in the lobes is considered in Sec 3, with applica- 
tion to high-power jets in Sec 4. The evolution of low- 
luminosity radio sources is discussed in Sec 5. 

2. GENERIC MODEL FOR JET-FED RADIO SOURCES 

We consider a generic model of extragalactic ra- 
dio sources, in which a relativistic jet emanates from 
the central engine, inflating a large diffuse radio 
emission region. In the current models of FR lis 
(jBlundell. Rawlings fc Willottl fl999h . it is generally as- 
sumed that particles are accelerated by shocks. The ac- 
celerated particles diffuse through the surrounding re- 
gion, referred to as the head region where hot spots are 
located, and are injected into the lobe. One assumes that 



these basic ingredients are also true for FR Is (cf. Sec 
5) except that for low-power jets the shocks occur close 
to the nucleus and the jets beyond the shocks continue 
to decelerate. A schematic diagram for a low-power jet 
is shown in figure [U The region where shocks occur is 
referred to as the flare region. In observations, the bright 
spots in the flare region can be regarded as the hot spot 
counterparts when compared to high-luminosity sources. 
The diffuse emission region, which extends well beyond 
the shock, is regarded as the lobe in analogy with FR lis. 

2.1. Self-similar expansion of jets 

For radio sources driven by high-powered jets, the 
effect of the pressure of the surrounding medium on 
the expansion of the lobes is negligible. As a re- 
sult, the jet-driven lobe expands in a self-similar man- 
ner. The distance from the central engine to the 
shock can be derived using the self-similar expan- 
sion argument (IFedorenko fc Zentsovalll986t IFalld I199TI: 
IKomissarov fc FalleHl998fh It is convenient to define a 
working surface at the shock, separating the jet from 
the lobe. Such hypothesis may appear to be oversim- 
plified, especially for multiple shocks that may occur, 
but it is useful in formulating a global model that links 
the various components of the jet-lobe system based en- 
ergy conservation. The input power by the jet across 
the surface is Qj and remains constant over the lifetime, 
where Qj is the jet power. The work done by the jet 
against the ram pressure is W ~ p9 2 v 2 r^, where r s is 
the radial distance to the surface, v ~ r s is the advanc- 
ing speed of the jet and is the opening angle of the 
jet. The density profile of the ambient medium, i.e. the 
halo of its parent galaxy, can be modeled as p — p c (l + 
r 2 /r 2 ) _/3 / 2 , where p c ~ 1.7 x 10~ 23 kgm~ 3 is the core 
dens ity and r c is the core r adius . typically about a few 
kpc (jGarrington fc Conwavlll99lHMulchaev fc Zabludofll 
119981: IBlundell. Rawlings fc Willottl I1999D . The profile 
index /3 is between 1 and 2.5. Writing W ~ Qjt, one 
obtains r, as a function of time: 



p c r'c 



1/(5-/3) 



£3/(5-/3) 



(1) 



with ci = [(5 - /3)/26»] 2/(5 - /3) and r > r c . For r < r c , 
one has r s (t) cx i 3 / 5 . 

For high-power jets, applicable for FR II sources, 2r s 
can be identified as the linear size of the lobe, denoted 
here by D{t) = 2r s (t) with r s (t) given by ((T|). One ob- 
tains 

/ x 3/(5-/3) 

D(t) = D l-) , (2) 

where D = 2c 1 (Q j /p c r^ )V(5-/3)f3/( 5 -« and t is the 
initial time when the size is Do- The initial size Dq is 
usually set to the radius of the hot-spot region. 

2.2. Pressure-limiting expansion 

The self-similar model considered in Sec 2.1 may not be 
applicable for low-power jets. The pressure in the lobes 
inflated by the low-power jets can drop rapidly to about 
the pressure of the external medium. Since the density 
decreases from the core r > r Cl the pressure of the ex- 
ternal medium decreases. We assume that the pressure 
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Fig. 1. — A schematic diagram for a radio source fed by a low- 
power jet. The flare region has an opening angle 0l much larger 
than the jet opening angle 9q. The dashed line indicates the outer 
boundary of the diffuse emission region that separates it from the 
external medium, i.e. the galaxy halo or intergalactic medium 
(IGM). In contrast to high-power jets, no shock develops at the 
boundary. 



in the lobes continues to balance the external pressure 
that decreases outward. We refer this scenario to as 
the pressure-limiting expansion. The argument for the 
pressure- limiting expansion can be understood as follows. 
Assume that the lobe expansion starts with near pres- 
sure balance pi ~ p cx . If a significant imbalance pi > p cx 
develops as a result of substantial drop in the external 
pressure as the lobe expands into the underdense region, 
the expansion will accelerate and the system will relax to 
near pressure balance on a time much shorter than the 
age of the source. One can derive the r elevant expan- 
sion law from the energy conservation (Bickncll 1986; 
Eilck & Shore 1989). Denoting the total energy of the 
lobe by £), the energy equation is given by 

dV 
dt' 



dt Qj 



Pi- 



tt) 



where one ignores the radiative losses and assumes that 
the injection power is Qj. It should be noted that pi 
is sum of all components (radiating particles, nonradiat- 
ing particles, magnetic field). Eq ([3]) implies that the 
main cause for decrease in the internal energy in the 
lobe is the volume work done against the external pres- 
sure. Assuming that the temperature of the external 
medium is constant, the external pressure can be writ- 
ten as p cx ~ p c (r c /r)' 3 . One may write the external 
pressure at the core as p c = uqUbTq ~ 1.4 x 10~ n Pa for 
= pc/mp — 10 5 m~ 3 and To = 10 7 K. As we are inter- 
ested in the pressure-limiting expansion with pi ~ p cx , 
one can substitute pi ~ Pcijc/rf f° r ® to obtain 



with 



D n = 2r r 



/r-n 


f 3 " 




\rlPcXJ 


V3r 





1 1/(3-/3) 



Qjt 



(4) 



(5) 



In deriving (@|, we assume Ei = piV/(T— 1) and V — x 7 " 3 , 
where T is the adiabatic index of the plasma in the lobe. 
The geometry of the lobe is described by the parameter 
X- Since there is no strong observational evidence for the 
dependence of the aspect ratio on the size, it is justified 
to treat x as an independent parameter. As x — 47r/3 
corresponds to a sphere, one generally has x ^ 47r/3. 

2.3. Constant-pressure expansion 

A special case is t he constant-pr e ssure expansion 
(jEilek fc Shore) ITTjsl IKaiser fc BesH I2007T I. With 
p ex =const, one obtains 



D(t) = D a 



1/3 



Dr, = r r 



in r-1 



3r*PexX r 



1/3 



(6) 

Eq (J5]) implies that the size of the radio sources driven by 
low-power jets grows much more slowly than the sources 
driven by high-power jets. 

2.4. Radio powers of lobes 

The total radio power as a function of time can be writ- 
ten as an integration of the single-particle power, P s , over 
the spatially-integrated, time-dependent particle distri- 
bution, iV(7,t), where 7 is the particle's Lorentz factor. 
In practice, one may assumes that each particle emits at 
the characteristic frequency v c — (3/4) vbJ 2 , producing 
a power spectrum P s 5{v — v c ), where vb is the gyrofre- 
quency and P s is the single-particle synchrotron power 
averaged on the particle's pitch angle, given by 

P s w I<7 T cUb1 2 , ( 7 ) 

where ctt is the Thomson cross section and Ub is the 
magnetic energy density. If the pitch angle distribution 
is maintained in the isotropic state throughout the evo- 
lution, the total spectral power can be expressed as 



Pu(t) = J P a 8(y-u c )N(rf,t)chf 



(8) 



where 7* = (4i//3^s) 1 / 2 is the Lorentz factor of parti- 
cles that emit synchrotron radiation predominantly at 
frequency v. Eq ([2]) can be used to eliminate £ in (JSJ) , 
giving rise to P v as a function of D. All the variables in 
([5]) are global, in the sense that they can be regarded as 
spatial averages. 

One can show that the synchrotron output (in radio) 
only comprises a tiny fraction of the total power input 
by the jet. The total number of relativistic particles that 
emit synchrotron radiation at frequency v is ~ nQjt"/* p , 
where 77 < 1 is the fraction of the jet power into the rela- 
tivistic particles, p is the particle spectral index and t is 
the age of the source. Since the total synchrotron power 
is ~ P v v, the ratio of the synchrotron power to the jet in- 
put (rjQj) is estimated as P v v/rjQj ~ o tC u B t(y j VB ) X I' 2 . 
For B = 50 nT, t = 1 Myr, and v = 1 GHz, this ratio is 
~ 10~ 4 , which implies that only a tiny fraction of the jet 
power is converted to synchrotron radiation. 

2.5. Limiting flux 

The instrument sensitivity places a lower-limit on the 
observable flux density. As a result, the P v -D tracks 
have a cut-off at which the flux density is too low to be 
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observable. To estimate the cut-off, one may write the 
flux density as 



F v ~ P v j 1 S > f„*, 



(9) 



with F v * the limiting flux density and S the effective 
surface area of the lobe. Hence, the following expression 



(10) 



defines a cut-off line. The sources below the cut-off line 
are undetectable due to the limit of the sensitivity. The 
simplest case is S ~ D 2 , which yields P v ~ F v1e D 2 . 

3. EVOLUTION OF THE EMITTING PLASMAS 

The temporal evolution of the relativistic particle spec- 
trum due to both adiabatic and radiative losses can be 
derived b y the usual me thod, i.e. solving a diffusion-loss 
equation (|Longairlll994fl . Since there is no strong evi- 
dence for frequency dependence of the source size, one 
may ignore spatial diffusion in the treatment of the evo- 
lution of the particle spectrum and only consider the 
spatially integrated distribution N(^,t). Here we out- 
line the main res ults based on the solutions discussed in 
iKardashevl (|1962h and concentrate on the case of the syn- 
chrotron losses in magnetic fields that decay with time. 

3.1. Time- dependent particle spectra 

When spatial diffusion is neglected, the diffusion-loss 
equation is simplified to the usual continuity equation, 

dN(j,t) d 



dt 



®1 



[bN( 7 ,t)]=q l ( 1 ,t), 



(11) 



where ©(7, t) is the particle injection rate, 6(7) = —7 is 
the energy loss rate, given by 



/,i,-| = (S^ + l),. 



_ 's'ics 



(12) 



where r s = 3m e c 2 / ' {AutcUb), t ics = 3m e c 2 /(4aTcU p h), 
C/ph is energy density of the seed photons. As usual, one 
assumes that the injection rate is time-independent, with 
a power-law distribution in energy: 

qi = <2o7~ P > 7i < 7 < 7m- (13) 

The injection rate is set to zero outside this range (i.e. 
qi = for 7 < 71 and 7 > 7 m ). For diffusive shock 
acceleration, the power index is close to p ~ 2. As- 
suming that the injection power is rjQj with r\ < 1 
the efficiency for the jet power going into the rela- 
tivistic particles, one obtains go = vQj/{(l) m e c2 ) = 

vQj(P - 2)/[m e c 2 (7i _p - 7 2 r p )], where (7) is the av- 
erage Lorentz factor of the relativistic particles. Two 
simplifications in regarding the continuity equation are 
made here. First, the effects of particle pitch angles can 
be eliminated by assuming that the pitch-angle distri- 
bution is isotropic and that the isotropic distribution 
can be maintained by efficient pitch angle scattering by 
plasma turbulence in the lobe. This assumption is valid 
if the isotropisation tim e is much shorter tha n the syn- 
chrotron cooling time (Daffe fc Perolal 119731 ) . Second, 
particles with 7 < 71 are ignored in the calculation of 
the power spectrum. Due to energy losses, particles 
can migrate from above the cut-off to below the cut- 
off (< 71). Thus, the particle spectrum that is initially 



TABLE 1 

EXPANSON LAWS 



Source 


Model 






ay 


FR I 


PLM 


0/(3 " 0) 


1/3 




FRII 


FEM 


(4 + /J)/(5-j8) 


(4 + 


/3)/r(5-/3) 



PLM— the limiting pressure model. FEM— the free-expansion 
model. The constant-pressure model can be regarded as a 
special case of PLM with /3 = 0. The parameters /3 and F 
arc given in Tabic [2] 

zero below the lower cut-off can become nonzero. We 
assume that the lower cut-off is the order magnitude of 
the bulk Lorentz factor of the jet, typicaly about ten. 
The typical Lorentz factor for high-frequency emission is 
7* - 10 3 (^/lGHz) 1 / 2 (50nT/B ) 1/2 - Since low- energy 
particles do not contribute to the high-frequency emis- 
sion, we can safely ignore the particles with 7 < "fi in 
our calculation. In practical calculations (cf. Sec. 4 and 
5), the results are not sensitive to the choice of j m as 
long as it is well above that required for the maximum 
observing frequency. 

Assuming that a particle is injected at to with the ini- 
tial Lorentz factor 70, the solution to the equation of the 
single particle's energy loss rate is obtained as 



7o 



l+7oV>(Mo)' 

' ^ dt' 



n 



(14) 
(15) 



where one assumes that the volume expands as V — 
Vo(t/to) 3av with Vq the initial volume at time to. The 
integration in (|15[) can be carried out provided that a spe- 
cific model for evolution of magnetic fields in the lobes 
is given. If magnetic fields are completely tangled, they 
can be treated as fluids. When the magnetic fields are in 
equipartition with particles, the density of magnetic en- 
ergy decreases with time, written as Ub = UBo(t/to)~ aB , 
where U bo is the initial density at to and the index as is 
a model-dependent constant. A summary of as and ay 
is given in Table [TJ The expansion law for magnetic fields 
is obtained by assuming equipartition with particles. In 
the limiting-pressure model, ag is obtained using @. 
Eq (fT5|) can be expressed as 



4>(t,t )=t 




, (16) 



(1 - a v )r IOS 

where r s o = 3m e c 2 / '(AutcU bo) and a 1 = as + ay. 

3.2. Characteristic ages 

The evolution of the emitting plasma in the lobes is 
characterized by three time scales: the adiabatic loss 
time, denoted by T a , the synchrotron time r s /7 and the 
ICS time t ics /j, where r s and t ics are defined in (TT2"j) . 
For ICS, one only considers the CMB radiation as the 
seed photons. The three time scales determine relative 
importance of these three energy loss processes in the 
evolution. 
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In the early phase both adiabatic and synchrotron 
losses can be important. Setting r a = r s /7, one obtains 
the characteristic age that separates the two energy loss 
regime, estimated as 

4/(4-3a B ) / 3 „ \ 2/(4-3a B ) 



ta 



tn 



to 



(17) 



where the Lorentz factor is assumed to be 7* = 

(4i//3i/bo) 1/2 (*Ao) q1b/4 . For a B > 4/3, the synchrotron 
losses dominate initially and at t > t a the adiabatic losses 
overtake the synchrotron losses. For as < 4/3, the adi- 
abtaic losses dominate first and then the synchrotron 
losses become important at t > t a . 

One may compare the ICS losses with the synchrotron 
losses. The energy loss rate due to ICS is j/t ics ~ 
^a T cU CMB /{3m e c 2 ), where U CMB = 3.6 x 1(T 14 (1 + 
z) 4 J m -3 is the CMB energy density at a redshift z. One 
may express U CMB in terms of an effective magnetic field 
Bcmb ~ 3.2 x lCT 10 (l + z) 2 T. Equating the ICS energy 
loss rate to the synchrotron loss rate yields a character- 
istic age: 



H — to 



B, 



CMB 



2/a B 



(18) 



The ICS losses become dominant over the synchrotron 
losses at t > t b . 

The characteristic age at which the adiabatic phase 
switches to the ICS phase is estimated to be 



t c 



to 



to 



4/(4+a B ) 



3^ \ 2/(4+Qf 



(19) 



Similar to the derivation of (fI7|) . the derivation of (|19|l 
involves replacing 7 by 7*. 

3.3. Analytical solutions 

We assume that particles are injected at a constant 
rate with spectrum qi — <7o7~ p and the initial condition 
N(p/,to) = 0. The formal solution for (|lip is written 
down in Appendix. There are two limits in which the 
exact analytical forms are well known. The first limit 
is when the adiabatic losses are dominant . One may set 
Ti — > 00, which leads to the exact solution (Eilek & Shore 
1989), 



N(rr,t) = 



<7o7 



-n 



1 + (p - 1)« 



1 * 



(p-l)a v +l 



.(20) 



Since the second term in the square brackets is generally 
much smaller than 1, Eq (|20p is approximately oc 7~ p , 
which implies that the spectral slope is not affected by 
the adiabatic losses and the whole spectrum raises pro- 
portionally with time. 

The second limit is when the energy losses are due to 
ICS of the CMB radiation or due to synchrotron losses in 
constant magnetic fields. One h as the well-known form 
(IKardashevlll96l IMelrosdflMI : 



N( 7 ,t): 



ToqoT 



p-1 



1-1- 



TQ 



p-1' 



> (21) 



where To = t ics for the ICS losses and tq — t s for 
the synchrotron losses. Eq ([21]) must be subject to the 



condition t — t < t /j. When 7 <C T~o/(t — to), i.e. 
the cooling time is much longer than the age, one has 
N(j,t) « qo(t — io)7~ p - This low energy limit can eas- 
ily be understood from the continuity equation in which 
the time derivative term, which describes the temporal 
evolution of the particle spectrum, is more important 
than the convection (in 7) term that corresponds to ra- 
diative cooling of the emitting particles. Thus, the par- 
ticle spectrum is cx t and its shape do not change. In 
the opposite limit in which the cooling time is much 
shorter than the age, 7 > To/(t — to), the convection 
term in (|11[) is more important than the time derivative 
term. One can directly integrate pip over 7 to obtain 
N(y,t) = <7o T o7~ p-1 /(p — !)• Since for radio sources, it 
is usually true that the radiative cooling time is consid- 
erably shorter than the source age, the particle spectrum 
in the high energy approximation is the more relevant. 

Apart from these two special cases, a third case of rel- 
evance, especially for low-luminosity sources (cf. Sec 5), 
is the synchrotron losses in magnetic fields that slowly 
decay with time. There is no simple analytical solution, 
though the partic ular case where ctg < 1, p = 2 and 
3 was discussed in [Eilek fc Shore (1989). However, one 
can express the formal solution in terms of the hypcrgco- 
metric function. The derivation is outlined in Appendix. 
Here we only discuss the case c*b = 1, which corresponds 
to (3 = 3/2 in the pressure-limiting expansion (cf. Sec. 
5). Other examples are considered in the Appendix. The 
solution for as — 1 has the form, 



N( % t): 



<iotj 



M(p-l, P ,l/0-(l-£M*/*b)r 



xM(p-l,p,l/£-ln(t/to)) 



(22) 



is 



where £ = 7to/ r sO < l/hi(f/fo) an< 4 M(a,b,x) 
the Kummer's function which has the asymptotic prop- 

erties M(p — l,p,x) ~ [p — l)x~ 1 e x for 1 > 1 

(|Abramowitz fc Stegunlll965l ). Thus, in the low energy 
limit £ ln(f/fo) <C 1, one has 



N(-/,t)^q tr 



1 - 



to 



7*0 



■In 

T S V to 



p-2 



, (23) 



which is similar to the low energy limit of (|2"Tj) . In the 
high-energy regime, one has (cf. Appendix) 

JV( 7s t)« ®Ztt-P-iM(p-l,p,l/0. (24) 
p-1 

The spectrum with an initial power-index p steepens to 

4. HIGH-LUMINOSITY RADIO GALAXIES 

Since FR lis have already been considered in the lit- 
erature, here we rederive the main features of their 
evolutionary tracks for nearby high luminosity sources 
(z <C 1), in particular the 'knee' feature due to the tran- 
sition from the expansion dominated by adiabatic losses 
to that by ICS of the CMB. To model the radiative evo- 
lution of the lobe we ignore the details of how relativistic 
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Fig. 2.— P„-D tracks for FR lis at 1.4 GHz 



top to bottom correspond to 



2 X 10 40 W 



The plots from 
5 X 10 39 W and 



10 39 W, respectively. The initial size and core radius are assumed 
to be 2.5 kpc. 



particles are injected into the lobe. One model for the 
injection is that the particles are accelerated in the shock 
at the hot spots. The accelerated particles diffuse across 
the head region where they are subject to synchrotron 
losses. Since the pressure in the head region remains ap- 
proximately constant as suggested from observations, the 
accelerated particles are subject to severe adiabatic losses 
as they enter the lobe whose pressure is a decreasin g func- 
tion of time (|Blundell. Rawlings fc Willottl fl999H . As a 
result, the radio power declines more rapidly than that 
inferred from observations. One possible remedy for this 
enhanced loss is that particle re-acc eleration is ongoing in 
the lobe (jManolakou fe Kirkl 120021). This would l ead to 
similar results to that obtained bv lKaiser fc Bestl (2007) 
based on the time-independent injection. Here without 
going into a specific injection model, we adopt the sim- 
ilar assumption that the injection is time independent, 
with a power-law energy distribution. 

The evolution of the emitting plasma in the lobes can 
be characterized by three separate phases: 1) the ini- 
tial build-up phase in which the total energy in the lobe 
increases with time due to the injection of both kinetic 
energy by relativistic particles and magnetic energy, 2) 
the adiabatic phase, and 3) the ICS phase. The rele- 
vant times t a , tb and t c can be estimated in a specific 
model f or the lobe pressure pi. We follow the proce- 
dure in iBlundell. Rawlings fe Willottl (|1999H by setting 
the pressure to that downstream of the bow shock, giv- 



ing pi oc t 



-(4+/3)/(5-/3) 



We also ignore the large-scale 



ordered magnetic fields and set the magnetic pressure 



to pi. The adiabatic expansion piV 



a v = (4 + 0)/T(5 - /3) (|Kaiser et alJH99l . One has 

/ . \ 33/5 

t a « 6.1 x 10- 3 ( ) Bf'\\^ Myr, (25) 



0.1 Myr 



t b « 62 



tc. 



0.1 Myr 



B 



14/11 
50 



(1 + *) 



-28/11 



Myr, 



(26) 



tc. 



135 



to 



S 6 36 Iyr 0.36 (1 + 2) -1.4 Myr; 



0.1 Myr , 

(27) 

where one assumes that B 50 = Bn/(50nT), v\ — 
f/(lGHz) and (3 = 3/2. That t b > t a implies that 
synchrotron losses are dominant only in the very early, 
build-up phase of the evolution. 

For high-power jets, by substituting (|19l) for one 
estimates the characteristic size at which the ICS losses 
dominate. This corresponds to a break or 'knee' at 



1% I '- 

to 



3/(5-/3) 



(28) 



Since in the adiabatic regime, the particle spectrum re- 
mains the same as the injection spectrum, the total spec- 
tral power for D < D c is P v oc Ug +P ^ 4 t, where the 
power-law index of the electron distribution, p, is as- 
sumed to be constant. Since t oc 

£)(5-/3)/3 ; the p v _ D 

relation is derived to be 



P„ oc D~ 



(5-/3)[(l+p)a B -4]. (29) 



Eq d29j reproduces (3) in iKaiser fc Bestl (|2007f ) when 
p = 2. One obtains 6 = (7/3 - 8)/12 « 5/24 for (3 = 3/2. 

When energy loss due to ICS is important, the number 
of particles with large Lorentz factors decreases leading 
to steepening of the particle distribution with the in- 
dex increasing to p + 1. The particle spectrum becomes 
time-independent, which is in contrast to the low-energy 
regime or the adiabatic regime in which the particle spec- 
trum increases with time. In the ICS regime (D > D c ), 



the total power is P v oc U 



(2+p)/4 
B 



D , with an index 



5= M5-p)(2+p)a B « 1. 



(30) 



const gives 



The i ndex is similar to that found in Kai ser fc Bestl 
(2007) (their Eq 7). The approximation is obtained when 
p = 2 and (3 = 3/2. As a result, the P v -D power-law in 
this regime is significantly steeper than (|29[) . 

Figure [2] shows three P v -D tracks at 1.4 GHz, with 
three different input powers. Here both the initial size 
and the core radius are taken to be 2 kpc. The tracks 
are obtained using the analytical solutions (f2"o) and (|2ip 
with 6 = 0.5, ry = 0.5, p c = 1.7 x lO^kgm" 3 , /3 = 3/2, 
and r = 4/3 (which is appropriate for relativistic plas- 
mas). We assume the particle spectrum to be the typ- 
ical one from the standard diffusive shock acceleration, 
characterised by a power- law with an index p = 2.1 and 
the lower- and upper-cutoff, (71,7™) = (5, 10 7 ) (cf. Eq 
[13]). The magnet ic field at the hot spots is assume d 
to be B = 30 nT (jBlundell. Rawlings fc Willottl [19991) . 
Since the initial particle spectrum is assumed to be zero, 
the spectral power, P„, increases rapidly to the phase 
when the spectral power reaches the maximum and the 
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adiabatic losses set in. During the initial phase, the syn- 
chrotron losses compete against replenishing of new par- 
ticles from the injection, with the latter dominating. The 
radio power decreases slowly as the source size grows, pri- 
marily due to energy losses through adiabatic expansion; 
during this phase, the power spectrum is P v ~ v~ a with 
a = (p — l)/2. When the size exceeds the characteristic 
size D c , the power starts to decrease rapidly due to the 
ICS losses. Since D c oc D% 28 with (3 = 3/2, the charac- 
teristic size D c is not particularly sensitive to the initial 
size Dq. As shown in the figure, the location D c of the 
break, i.e. the 'knee' feature remains roughly the same 
for different Dq- In the ICS regime, the power spectrum 
steepens to a = p/2. The tracks are the most strongly 
affected by the input power i]Qj, the magnetic field at 
the hot spot, -Bo, and the index for the density profile, 
p. The track height increases if one increases the input 
power or the magnetic field at the hot spot or both. In 
general, a larger (3 leads to steepening of the track slope 
in both the adiabatic regime and ICS regime. 

5. LOW-LUMINOSITY RADIO GALAXIES 

The pressure-limiting expansion model is applied to 
FR I sources. A low-power jet may undergo free expan- 
sion initially, but such initial phase can only last very 
briefly. Therefore, one can ignore this phase and assume 
that the jet expands slowly in the pressure-confined en- 
vironment. 

5.1. Observational data 

To test the models, we assembled a complete sample 
of low-luminosity radio galaxies from the northern zone 
of 2d F Galaxy Redshift Survey (2dFGRS) (jColless et alJ 
l2001f ) which has 1.4 G Hz radio continuum d ata available 
from both the NVSS (ICondon et al.l H998T) and FIRST 
(IBecker. White fc Helfandl 11995ft surveys. The 2dFGRS 
radio sample was selected by matching the NVSS and 
2dFGRS catalogue s using the techniques described by 
Sadl er et al.l (|2002f ). We then restricted our final sample 
to galaxies which satisfied the following criteria: 

• Classified as an AGN on the basis of the 2dFGRS 
optical spectrum (i.e. star-forming galaxies are ex- 
cluded) . 

• The total NVSS 1.4 GHz flux density is at least 
lOmJy, to allow an accurate measurement of the 
radio-source angular size. 

• The galaxy is also detected in the FIRST catalogue, 
which has higher angular resolution (with a 5 arcsec 
beam, compared to 45 arcsec for NVSS). 

5.2. Radio luminosity and source size 

This selection produced a final sample of 375 low- 
luminosity radio galaxies, with redshifts in the range 
z = 0.02 to 0.3 (median z = 0.136). All 375 sources 
were spatially resolved in the FIRST survey, so there are 
no upper limits on the angular size of these radio sources. 

We calculated the 1.4 GHz radio power and largest lin- 
ear size (LLS) of each source from the measured NVSS 
flux density and the largest angular size (LAS) mea- 
sured from either the NVSS catalogue or (if unresolved 
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Fig. 3. — Local low-luminosity radio galaxies at 1 .4 GHz, from 
the 2dFGRS sample described bv lSadler eFaTI J2003) . 



in NVSS) the FIRST catalogue. We adopted a cos- 
mology with Ho = 71 kms -1 Mpc" 1 , Sl m = 0.23 and 
£1a = 0.73 for these calculations. In the few cases where 
a source was resolved into two or more components in 
NVSS or FIRST, we measured the LAS across all compo- 
nents. The 1.4 GHz radio power of galaxies in this sam- 
ple spans a range from 10 22 to 10 26 5 WHz _1 (median 
10 24 - 2 WHz -1 ), and the LLS values range from 0.4kpc 
to 715 kpc (nedian 51 kpc). 

The radio P u -D diagram for the sample is shown in Fig- 
ure^ There are two interesting features that can be used 
to constrain the modeling. First, there are fewer sources 
with a small size in the high luminosity region. Second, 
there is a deficit in the number of sources with large sizes 
in the low luminosity region. The second feature may be 
understood qualitatively as due to the flux limiting ef- 
fect (cf. Sec. 2.5), i.e. because of the lower limit on the 
observable flux, the P v -D have a cutoff P v > F uie D 2 . 
The first feature suggests that the source, as we argue, 
expands at a rate considerably faster than the constant- 
pressure expansion (see Sec 5.2). As a result of increased 
losses, the initial increase in the luminosity is more grad- 
ual than that predicted by the constant-pressure expan- 
sion model. 

5.3. Power-size tracks 

We consider the scenario of constant-pressure expan- 
sion first. The external pressure p ex is held constant with 
Pi ~ Pex- We show that the constant-pressure model 
would generally underpredict large-size sources. The 
model overpredicts small-size sources with high luminosi- 
ties. If equipartition applies, one has pb ~ Pi ~ const. 
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TABLE 2 

Model parameters for low-luminosity sources 



(irtm nnn tint 




T)o~n ti iti run 




Jet 


Qj 


Injection power (W) 


10 32 - 5 x 10 36 




V 


Efficiency of power injection 


0.5 


Lobe 


X 


Geometric factor 


I 




r 


Adiabatic index 


4/3 




p 


Relativistic particle spectral index 


2.1 




71 


Minimum Lorcntz factor 


5 




7m 


Maximum Lorentz factor 


10 7 


Environment 


Pc 


Core density (kgm -3 ) 


1.7 x icr 23 




r c 



Core radius (kpc) 

Power index of the density radial profile 


1.5 - 3 

1.5 (PLM), (CPM*) 




To 


Temperature (K) 


10 7 




Pc 


Core pressure (Pa) 


3 x icr 11 



CPM— the constant pressure model. Since p c 
independent parameter in our model. 
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Fig. 4. — P v -D tracks for low-luminosity radio galaxies. The 
tracks are obtained in the constant-pressure expansion model. 
Plots from top to bottom correspond to (Qj, Do) = (5 X 
10 36 W, 3 kpc), (10 36 W, 2 kpc), (10 35 W, 1.5 kpc), respectively The 
tracks increase much more slowly in the initial phase compared 
with FR lis. The plots are terminated at t = 10 4 Myr. 



Setting ay = 1/3 and as = 0, one has 

' q \ 1/2 
3^B0 \ 

4r/ / 



ta — oT s Q 



0.4R 



3/2 -1/2 



v\ ' Myr. (31) 



Since tb — > oo, the synchrotron losses are always domi- 
nant over the ICS losses. The synchrotron losses domi- 
nate over the adiabatic losses at t > t a . Eq (|3"Tj) implies 
the characteristic size 



T S 



l.6D B 



-1/2 -1/6 
10 "l 



1/6 



tr. 



0.1 Myr 



-1/3 



(32) 



— Pc^bTq /nip, the pressure p c is treated as an 

One may express the total power as a power-law of the 
size, P v ~ D~ s . For D < D a , one has 5 = —3. Since 
the particle spectrum in the synchrotron regime with 
t <C to + t s o/7* has the same slope as in the adiabatic 
case, one has 5 = — 3 for D > D a as well. The parti- 
cle spectrum due to the synchrotron losses in constant 
magnetic fields has the same form as (|2~T1) with tics re- 
placed by t s q. The power grows with increasing size. 
However, when t ~ to +t s o/7*, corresponding to the size 

~ D [l + (r s0 Ao)(3^Bo/4^) 1/2 ] 1/2 , the spectrum of the 
emitting particles steepens from p to p + 1 and becomes 
independent of time. As a result, the power remains 
constant with S = 0. The model would predict a steep 
spectrum with a = p/2 = 1 for p — 2. Although there 
is no direct measurements of the radio spectral index 
distrbution in the sample of 2dFGRS galaxies used to 
derive the P v -D diagram in Figure [3j there is strong evi- 
dence from the lower frequency data that low-luminosity 
radio galalxies on average have flatter radio spectr a 
(jMauch et al.l[2fM [Bhindell. Rawlings fc WillottlllQQflh 
A possible interpretation for this spectral trend is that 
for low-luminosity radio galaxies, t he core component be- 
comes important (|Slee et al.lll994D . i.e. one may obtain 
flat spectra by appealing to self-absorption. However, the 
steep spectral feature predicted by the constant-pressure 
model would extend to the majority of low-luminosity 
sources, which is not consistent with observations. 

Figure S] shows plots of P v as a function of D for low- 
luminosity sources (FR Is), obtained in the constant- 
pressure model with z = and (3 = 0. The model param- 
eters and their typical values are listed in Table [21 If one 
considers the external medium as an ideal gas, the three 
parameters, p c , Tq and p c are not indepedent. Here we 
treat the external pressure as an independent parameter 
and is assumed to be p ex = p c — 3 x 10~ n Pa. As we 
assume that the magnetic energy is in equipartition with 
kinetic energy of particles, the magnetic field pressure 
is not an independent parameter. The justifications for 
the choice of the particle spectrum are the same as that 
for high-luminosity sources-particles are accelerated by 
diffusive shocks (cf. Sec. 4). The typical value of the 
adiabatic index, T = 4/3, is applicable for relativistic 
plasmas. (For a cold plasma, one has T = 5/3.) The 
dividing line between high- and low-luminosity sources 
in terms of the input power from the jet is ~ 10 37 W 
(|Ledlow fc Owedll996l ). The plots are obtained using 
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Fig. 5. 
0.02 Myr, 



20 40 60 80 100 

B(nJ) 

— Plots of t a (solid), t\y (dotted), t c (dashed) with to 
2 = and 8 = 3/2. 



three different input powers that are below this charac- 
teristic power. The dots represent the observations (cf. 
Sec. 5.2). The plots are terminated at t = 10 4 Myr. The 
luminosities remain constant after reaching their max- 
ima. The luminosities may eventually drop below de- 
tection. For example, one may invoke a model in which 
the lobes are assumed to expand into a much lower den- 
sity region, which can result in a rapid d ecrease in the 
luminosity (Copal-Kris hna fc Wiitalll988l) . However, as 
shown in figure 01 the constant-pressure model underpre- 
dicts the large-size sources while predicts excess number 
of the small-size sources. 
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Fig. 6. — As in figure [4] but obtained in the pressure-limitin, 
model with Q, = 1.5 X lo" W and 8 = 3/2. The solid and dashcc 



8 = 3/2. 



dashed 



lines correspond to r c = 2 kpc and r c = 1.5 kpc, respectively. 



5.4. Pressure-limiting approximation 

When the lobe expansion is pressure limited with j3 ^ 
0, the three characteristic times discussed in Sec 3.2 are 
now given by 



L«2x 10" 



t b « 2.4 x 10 d 



t c w 2.6 x 10 2 



tr 



0.1 Myr 



SsoV'Myr, 



(33) 



in 



0.1 Myr 



B 5 2 (l + *)- 4 Myr, (34) 



to 



0.1 Myr 



1/5 



B^ 7 i/r 3/T (l + *)~ 1 - 6 Myr > 



(35) 

where /3 = 3/2 and olb — PI (3-/3). Here the estimate for 
as is obtained assuming energy equipartition between 
the magnetic field and particles (cf. Tabled]). Plots of 
(t a , %,t c ) are shown in figure [5] Like the high-luminosity 
sources, both t a and tb cx Bq are strongly dependent 
on the initial magnetic field which is assumed to be the 
magnetic field in the flare region. At Bq > 26 nT, the 
energy loss processes dominate along the evolutionary 
track in the following order: adiabatic, synchrotron, ICS; 
in the low field case Bo < 26 nT, the ICS dominance 
occurs before the synchroron loss. It is interesting to note 
that tb is very sensitive to z compared to high-luminosity 
sources. For local sources (z -C 1), tb is over 100 Myr. 
Thus, for local sources, the ICS losses play a role only in 
the very late stage of the evolution. 

The Py-D tracks can be obtained analytically. There 
are three relevant phases: the initial rapid rise, followed 



P v (W/Hz) 



25 



10 



24 



10 



23 




10 



10 



10 



D (kpc) 



Fig. 7. — As in figure[6]for two different jet powers, Qj = 10 36 W 
(solid line) and 1.5 X 10 36 W (dashed line). We assume Dq = 2 kpc. 
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Fig. 8. — High-luminosity sources (dashed) and low-luminosity 
sources (solid). The dots (from left to right) on each curve indicate 
ages; for high-luminosity sources, they represent 0.1 Myr, lMyr, 
10 Myr and 10 2 Myr. For low-luminosity sources, the dots represent 
1 Myr, 10 Myr, 10 2 Myr, and 10 3 Myr. 

by the more gradual increase in luminosity, and then 
decline in luminosity. In the gradual increasing phase, 
the particle spectrum is AT (7, t) oc t^~ p . The P v -D track 
can be approximated by a power-law, P v ~ D~ s . From 
© using (jU), one finds 8 = -(3- /3)[4 - a B (l +p)}/4 = 
—3/8 for p = 2 and f3 — 3/2. When the self- absorption 
effect is neglected, the corresponding radio spectrum is 
a = (p — l)/2. When t ~ t a , the source enters the 
declining phase in which the magnetic fields in the lobe 
become too low and given a fixed observation frequency, 
the Lorentz factor of the emitting particles shifts to a 
much higher value. The time t a when the track starts to 
turn over is sensitive to the initial time to and magnetic 
field B . For to = 10 4 yr and B — 50 nT, one estimates 
t a 1=3 20 Myr. For t > t a , the synchrotron losses become 
dominant and the particle spectrum is stationary with 
the spectral slope steepening to p + A with A a function 
of f3. This leads to S = [p + A)/3/2. The radio spectrum 
steepens to a — (p + A — l)/2. For (3 = 3/2, one has 
A~l. 

Figure [6] shows the P v -D tracks in the pressure- 
limiting case with the parameters given in Table [21 Al- 
though there is a spread in j3, here we take the typ- 
ical value (3 = 3/2. The pressure at the core is as- 
sumed to be p c — 3 x 10~ n Pa corresponding the den- 
sity n = 2 x 10 5 m -3 and the temperature T = 10 7 K. 
As in figure 21 the plots are overlaid on the observa- 
tional data of local luminosity sources represented by 
dots. Figure [7] shows how the turn-over depends on the 
input power Qj. For low Qj, the track turns over at a 
much smaller size. The evolution pace is characterized 
by dD/dt ~ (Do/t )(t/t )- {2 - f:i)/{3 ~ f:l) . Since the initial 



size Dq is treated as an input parameter in the plots, 
([5]) implies that a small Qj corresponds to a long to, 
i.e. the source moves more slowly along the track than 
otherwise. Thus, the input power determines the evolu- 
tion pace along the track. As in the constant pressure 
case, we adopt equipartition in our calculation; Thus, 
there is a direct connection between the magnetic pres- 
sure in the lobe and the external pressure. Since the syn- 
chrotron power is oc Ug 1 ^' i , the track height depends 
more strongly on p c than the jet input power (Qj). One 
can calculate P V (D) for an intermediate case < f3 < 3/2 
as well. It can be shown that the resultant track increases 
rapidly in the initial phase which similar to the constant- 
pressure case but the radio power drops off rapidly at 
small sizes. This early steep drop off is due to that the 
particle spectrum become very steep with p+ 1/(1 — as) 
(cf. the comments below Eq [A10]). In the case of a 
steep density profile (j3 > 12/7), one expects the track 
to change to the shape with a declining trend. This can 
be understood qualitatively that the lobe expands more 
rapidly in an external medium with rapid decreasing den- 
sity in radial direction. 

Figure [5] shows comparison of FR Is (solid lines) and 
FR lis (dashed line). One assumes Qj = 10 40 W and 
(Do, r c ) = (2 kpc, 2 kpc) in plotting the dashed curve and 
Q-j = 10 36 W and (D ,r c ) = (1.5 kpc, 1.5 kpc) in plot- 
ting the upper solid curve and (Do,r c ) = (2 kpc, 2 kpc) 
in the lower solid curve. Except for the initial brief 
increase, the luminosities of FR lis decay during their 
life time. By contrast, our model predicts that low- 
luminosity sources grow with increasing luminosities. 
The growth of radio power stops when the synchrotron 
or ICS losses become dominant over the adiabatic losses. 
For Qj = 1.5 x 10 36 W, the radio power peaks at about 
100 Myr. 

6. CONCLUSIONS AND DISCUSSION 

We consider the evolution of low-luminosity radio 
galaxies with z -C 1. The radio power as a function of the 
source's size is derived based on a generic model that uses 
global parameters of the jet-lobe system. In the global 
model discussed here, we assume that the lobe inflated 
by a low-power jet undergoes pressure-limiting expan- 
sion, in which the system is approximately in pressure- 
balance against the external medium. Thus, the energy 
equation (Eq [5]) for the lobe can be solved by replac- 
ing the lobe pressure, pi, with the external pressure, p ex . 
The source's size as a function of time can be derived in 
both cases of constant external pressure and decreasing 
external pressure. In the calculation of the P v -D track 
we use the standard theory of nonstationary spectra of 
relativistic particles that are subject to both adiabatic 
and radiative losses. The predicted P V (D) can be com- 
pared with the P v -D diagram from observations. The 
main conclusions are summarised as follows. 

1) The pressure-limiting expansion model (with /? ~ 3/2) 
predicts P v -D tracks that are generally consistent with 
observations; in particular, it predicts fewer small-size 
sources in the high luminosity region in the radio P v -D 
diagram. Since the emitting particles suffer expansion 
losses, the luminosity increases gradually untill the syn- 
chrotron or ICS losses become dominant. 

2) The constant-pressure expansion model (with j3 = 0) 
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generally predicts an overabundance of small-size FR I 
radio sources with a relatively high radio power and at 
the same time underpredicts the large-size FR I sources. 
These features are not consistent with observations (cf. 
Figure HJ. 

3) In the pressure-limiting expansion model, the physi- 
cal conditions of the external medium, e.g. p cx , play an 
important role in the evolution. By assuming equiparti- 
tion, one expects to see a direct link between the radio 
power and the external pressure, i.e. the higher external 
pressure the higher radio power. The P y -D track would 
also turn over at a much smaller size. 

It is interesting to compare the model prediction 
for low-luminosity sources with that of high-luminosity 
sources. For low-power jets, the lobes expand slowly 
in near pressure equilibrium with the external pressure. 
Thus, FR Is evolve much more slowly compared with 
FR lis. Our model predicts that the radio power of FR 
Is may increase throughout most of their life time. By 
contrast, the radio power of FR lis declines in most of 
their life time (except for the brief initial rise) and can 
be well described by two power-laws. The initial rise for 
FR lis is very brief, less than 0.1 Myr for the parameters 
adopted in Figure EJ while for FR Is, this can last more 
than 100 Myr. These differences suggest that these two 
types of source evolve differently. 

In our global model, both the details of the particle in- 
jection and the effect of spatial diffusion are not treated. 
The effect of spatial diffusion can be important in par- 
ticle transport-this is particularly the case when parti- 
cle acceleration is confi ned to a localized site, say th e 
flare region ('hot spot') (|Eilek, Melrose fc Walkerlll997| ). 
The accelerated particles need to diffuse across the re- 
gion surrounding the acceleration site. Note that this 
requirement may be relaxed if particles are accelerated 
by multiple weak shocks distributed over an extended 
region in the lobe. Particle diffusion should strongly de- 
pend on plasma turbulence in the region concerned. To 



obtain the solution for particle diffusion, one needs to 
deal with the full diffusion-loss equati on that includ es 
particle diffusion in plasma turbulence (|Longairlll994l ). 

One should emphasise that the pressure in the energy 
equation © is the total pressure that may consist of ra- 
diating and nonradiating particles as well as magnetic 
pressure. For convenience, we consider only one species 
of particles (radiating particles). X-ray observations sug- 
gest that in some of the known FR I sources, the equipar- 
tition pressure in t he lobe is substanti ally lower than the 
external pressure ([Croston et all 120081 ) . One possible ex- 
planation is that the missing pressure may be provided 
by nonradiating particles due to entrainment. Although 
in principle, the case can be treated in a similar way 
to that presented here for single species of particles by 
choosing a low rj and T —> 5/3, An extension of this 
model to including multi-component plasmas in the lobe 
will be considered in future work. 

Finally we comment on the equipartition assumption 
adopted in the calculation of the radio power. The mag- 
netic fields are assumed to be completely entangled (i.e. 
we ignores the effect of the mean field). The equiparti- 
tion assumption is reasonable as recent Chandra X-ray 
observations suggest that the magnetic fields in lob es are 
close to the equipartition field (|Birzan. et alj |2008) . It is 
worth noting that the frozen flux argument would sug- 
gest that the magnetic energy density in an expanding 
lobe decreases with time much faster than cx l/t 13 ^ 3 ^^. 
To maintain equipartition, one requires eithe r that the 
magn etic turbulence be generated in the lobe (|De Youngl 
1980) or that the magnet ic energy be injected into the 
lobe (|Eilek fc Shore! fl989l) . 
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APPENDIX 

SYNCHROTRON LOSSES IN DECAYING MAGNETIC FIELDS 

Similar to the method discussed in lEilek fc Shore! ([19895 > one ma 7 solve (fTTj) by solving the ordinary differential 
equation: 

dN( l0 ,t) 



dt 



6 7 (7 ,£)N(7o,t) + qi{lo,t), 



(Al) 



where 70 is treated a parameter and the functions 6 7 (7o,i), N("f ,t) and 31(70, £) correspond to db/dj, N{y,t) and 
Qiil) with 7 written as a function of (70, t): 



7o 



1 + 7oX 



to 



t s0 (1 - a B ) 



1 — Of 



(A2) 



The solution is 



Co 



1 — OLB 



I(7o,0, /(7b, t) 



-> . SO ( i\—a a 

1 + \X - 1 

1 — OLB 



p-2 



dx, (A3) 



where £0 = 70^0/^0 and x = t/to- The final form N(j,t) can be obtained by transforming (|A3|) back to the domain 
(7,i), using 



7o = 7 



1 - OL B 



l-a B 



(x x - aB - l) 



(A4) 
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Special cases 

Consider two special cases a B — 1/2 and a B = 1- For a B = 1, one has 
iV(7o, h) = qotolo P (1 + £o In a;) 2 I(lo, t), 



(A5) 



(1 + Co Inxf- 1 M(p - l,p, (1 + £o lna:)/£ ) 



&(P-1) 

-M(p-l,p,l/e ) 

where M(a, b, z) is the Rummer's function, defined as (|Abramowitz fc Stegunl figGS') 



M(a, b, z) 



T(b) 



T(b - a)T(a) Jo 



oZ t t {a-l)^ 1 _ i y>-a-l dt 



(A6) 



(A7) 



Eq ([22]) can be obtained by substituting (|A4I) back to (IA5|) and (|A6|) . The low energy limit (|23|) can be obtained by 
tak ing the limit £ln(i/i ) <C 1. The high energy limit (l2~i|) can be derived by taking either the limit £ ln(i/*o) 00 
in (|A3|) or the limit £ ln(t/to) — > 1 in the final form. 
For a B = 1/2, the integral can be evaluated exactly as 



/(To,*) 



2£ Mp-1) 



This gives 



N( 7 ,t) 



qotoi p 



1 + 2£ (a a/2 - 1) 



2£pv^ - 1 - 



- p-l 






2£o ( 



2£o(1 + Cp-1)z 1/2 ) -1 



- (2CoP - 1) 



2p( P -i)e 2 v 

In the low energy limit ^(^/x — 1) <gc 1, one has 



1-2^-1) 



2^(p-l + Vi)-l 



1/2' 



-1 2 



(A8) 



(A9) 



(A10) 



One can obtain the high energy limit by setting 2£(y / a: — 1) — > 1 in (|A9|) . which leads to N( j ,t) s a 
(go/2p)r s o(T s o/io)7 _p_2 - The power index obtained here is consistent with that derived in lEilek fc Shore! fl989). 
It is worth noting that their approximation gives an index ~ p + 1/(1 — which is not applicable for = 1. 

General case 

Generally, the integral can be expressed in terms of the hypergeometric function using the following integral 



{l + At) a {l + Btydt-- 



(B AY 



o 



(1 + B) 1+0 F(l + /?, -a; 2 + /3; A 
A 



(1 + (3)B 1 + 

-F(l + ^-a;2 + / 3;-^) 
where the hypergeometric function is defined as (jAbramowitz fc StegunllT965j ) 

r(c) 



l + B 
A-B 



F(a, b\ c; z) 



t b - 1 {i-ty-'>- l {i- zt)- a db. 



z-b-l 



One obtains 



r(6)r(c-6) Jo 

[1 - (1 - a B )/to] 



(All) 



(A12) 



Co (P-l) 



(l + B)P- 1 J F(j»-l > -a;p; 



AQ. + B) 
A-B 



F(p-1, 



A 



-a;p; 



A-B 



(A13) 
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where 



A = x 1 ~ aB -l, B = 



1 - a B 



1) 



as 



1 - a B 



(A14) 



Again here th e fin al solution can be obtained by substituting (IA13|) into (|A3[) and transforming it back to the (7, t) 
domain using (|A4I) . This gives 



N( 1 ,t) = q Q t Ql -P 



1 - OLB 



p-2 



Alternatively, the solution of the continuity equation (fTT|) can be written as ()Kardashevl 1 1 9621 ) 



to 



IP' 2 ^\ 



with 



/ t ' v*'7 ri(t") : 
where T\(t) is defined in (JT5J). Eq (|A17j) can be integrated to yield 



1> = ti 



1 



(a' - l)r s0 V^o 



1 



1 



(ay - 1)t/cs 



The solution satisfies the initial condition 2V(7, 0) = 0. For «b = 1, one has 



^, t ') = -^lnfi 



(A15) 



(A16) 



(A17) 



(A18) 



(A19) 



Substituting Eq (1A"T9|) into (|AT6|) leads to 
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